Gordon Metric Revisited 
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We show that Gordon metric belongs to a larger class of geometries, which are responsible to 
describe the paths of accelerated bodies in moving dielectrics as geodesies in a metric q M „ different 
from the background one. This map depends only on the background metric and on the motion of 
the bodies under consideration. As a consequence, this method describes a more general property 
that concerns the elimination of any kind of force acting on bodies by a suitable change of the 
substratum metric. 

PACS numbers: 04.20.-q 



I. INTRODUCTION 

In 1923 Gordon [l| made a seminal suggestion to treat 
the propagation of electromagnetic waves in a moving 
dielectric, modifying the metric structure of the back- 
ground. He showed that the electromagnetic waves prop- 
agate as geodesies not in the background geometry rj^ 
but instead in the effective metric 
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rf v + {en - 1) v u , 



(1) 



where e and fi are constant parameters that characterize 
the dielectric and i> M is the four-velocity of the material 
under consideration (which is not necessarily constant). 
Later, it was recognized that this interpretation could be 
used to describe nonlinear structures even when e and fi 
depend on the intensity of the electromagnetic field |2J or 
more complicated functions of the field strengths [3j. In 
all these cases, the causal cone, which is associated to the 
effective metric, does not coincide with the null-cone of 
the theory. The origin of this modification is due to the 
presence of a moving dielectric, which changes the paths 
of the electromagnetic waves inside this medium. 

We then face the question: could such particular de- 
scription of the electromagnetic waves in moving di- 
electrics be generalized for other cases, in which acceler- 
ated paths due to other kind of forces would be described 
as geodesic motions in an associated metric? We shall sec 
that the answer is affirmative and this kinematical map 
depends only upon the acceleration of the body and the 
background metric. 

This method allows us to geometrize any force in the 
sense that arbitrary accelerated body in a given met- 
ric substratum g^ is equivalently described as geodesic 
motion in an effective geometry q^ v . We start by an- 
alyzing the generalization of Gordon metric concerning 
the propagation of electromagnetic waves inside arbitrary 
dielectric media. This procedure mimics the trajectory 



which led to the geomctrization of the gravitational field 
as it was done by general relativity (GR). This means 
to describe the effects of acceleration of a particle on a 
gravitational field by a change of the space-time met- 
ric, according to Einstein's approach. We compare this 
effective geometry to the metric in the post-Newtonian 
approximation in order to sec if it is possible to repro- 
duce some results of GR. Such procedure is exactly what 
happens in the analog models of gravitation that deals 
with systems kinematically equivalence, but dynamically 
distinct. 



II. GEOMETRIZING ACCELERATED PATHS 

Consider a vector field 1 with norm N = u a u a in a 
given background g^u . The acceleration of is given by 



Now let us construct an associated metric tensor q^ v , in 
which the vector satisfies the equation 



(2) 



where 1 1 means covariant derivative with respect to q^ v 
and f(p) is an arbitrary function of the parameter p along 
the curve. The contravariant components of the vector 



field are defined by m p 



and, consequently, the 



norm defined in is given hy N = q^ v u^u v . Whenever 
is cither a gradient or a normalized vector field, f(j>) 
can be set equal to zero without loss of generality. 
Developing Eq. (f2]), we obtain 



(3) 
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1 The natural velocity field we use to develop this section is a given 
1— form field u M . In Gordon approach, the wave vector fc M is a 



gradient and, therefore, an exact 1— form. 
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whore [] means skew-symmetrization. Choosing q lu , such 
that this equation is verified, then the accelerated path 
of u a in becomes a geodesic motion in the associated 
q^ v . Note that this is independent of its functional form. 
We will show that Gordon metric is a particular exam- 
ple of this procedure and that there is a class of metrics 
which play the same role in such "unforced-motion" pro- 
cess exhibiting different geometrical properties for each 
element of the class. 



III. LIGHT PATHS ON MOVING DIELECTRIC: 
GORDON APPROACH 

Let us define two skew-symmetric tensors and P^ v 
representing the electromagnetic field inside the material 
medium. These tensors are expressed in terms of the field 
strengths and H v and field excitations and B^ as 
follows 



P, 



where w M is a given four-vector comoving with the di- 
electric and rjfj, u a ^ is the Levi-Civita tensor. We assume 
that the electromagnetic properties of the medium are 
characterized by the constitutive relations 



D a — '(E,H)E V , B a — \ia '(E,H)H V , 

where e a "(E,H) and [i a u (E,H) are arbitrary tensors. 
Consider Maxwell equations on dielectric media Q with 
permittivity e and permeability \i that characterize the 
dielectric: 



= 0, 



"F^ ;v = 0. 



(4) 



Prom now on, we take the background metric as flat 
Minkowski space-time and assume that [i = /io is a con- 
stant and e = e(E), where E = y/-E a E a and E a is 
the electric field. It is straightforward to generalize these 
equations to arbitrary curved space-time. Indeed, sup- 
pose an observer with velocity comoving with the di- 
electric and such that v^-v = 0. Then, Eqs. (j4|) written 
in terms of the displacement vectors D M and B^ become 



D». v v v - D v . v v> i + if al3 v a H^ = 0, 
B»- V v v - B v - V v» - if vaP v a E^ v = 0. 



(5) 



The projection with respect to v M yields the four in- 
dependent non-linear equations of motion describing the 
electromagnetic field inside the dielectric medium: 



e E a - a — 6 E a .p — 0, 

HoH a ;a = 0, 

eE x - £&tPEt Ema + n^P-VpH^ = 0, 
,i ij A - v^VpE^p = 0. 



(0) 



We define the unitary vector by setting E^ = El^, 
where Z M satisfies l a l a = —1. 

We use Hadamard conditions [B| to obtain the prop- 
agation waves through the characteristics surface £ (for 
details, see appendix). The symbol [X] E represents the 
discontinuity of X through this surface. Then, the dis- 
continuities of Eqs. © become 



(7) 



where e M (x) and h^(x) are the amplitudes of the disconti- 
nuities and fc M = dfiY, is the wave vector. Thus, it follows 
that 



ek a e a - ^E a e a E' } kp = 0, 
[i h a k a = 0, 

ek a v a e>* - ^E x e x v a k a E^ + T)^ va ^k v V a hp = 0, 
fi k a v a h x - r} x ^ a k v p e a = 0, 



(8) 



where e' is the derivative of e with respect to E. Combin- 
ing these equations we obtain the following intermediary 
relation 

-^E x e x v a k a E^ = 0, 

(9) 

which multiplying by E„ yields the dispersion relation 



rT + (^oe - 1 + w'E)v»v v - j^E^E") k^k v = 0. 

(10) 

We see that the envelop of discontinuity propagates dif- 
ferently from Minkowski light-cone of the linear Maxwell 
theory. In this case, the causal structure is given by an 
effective Riemannian geometry 2 g^ v . From this point of 
view, k^ is null- like in g^ v ', namely, 



gTk^ky = o. 



(ii) 



2 Mathematically, the metric tensor is a covariant tensor of rank 
2. However, in this paper, we sometimes shall call "metric" a 
contravariant tensor of rank 2. In particular, that is the way the 
Gordon metric appears naturally. 
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The expression of the effective geometry is given by 



and 



tf» = if* + ( Moe - 1 + m > 'E)v»v v - —m v . (12) 
A simple calculation show that its inverse is 



9nu = Vfiu ~ 1 - 



where 



w ( lt!L ,V.t-V. (13) 



e 

In particular, when e is a constant, this formula re- 
duces to Gordon pioneer work, in which was shown that 
the waves propagate as geodesies not in the background 
geometry r\^ v but instead in the effective metric 



n 



(14) 



which depends only on the dielectric properties /io, e and 
The magnitude N of the wave vector in Minkowski 
space-time (written in terms of dielectric properties) is 
determined by Gordon relation 



g^knh = W + (cmo - 1) k M k v = 0, : 

N = (1- t* e){k.v) 2 , 



(15) 



where k.v = k a v a . 

The analysis of the wave propagation in material me- 
dia and the study of effective geometry are particularly 
interesting in the investigation of analog model 0, [f| for 
the understanding of kinematics! properties at very small 
scale of astrophysical objects (see details in @, Q ) • We 
quote Hawking radiation Q and Unruh's work on exper- 
imental black hole evaporation (Toj which are systemati- 
cally studied and the modeling of specific dielectric media 
is developed in order to eventually detect these tiny ef- 
fects. A more complete discussion on this topic was given 
in [ll|, [HI and references therein. Here we shall point the 
similarity between these geometries in a later section. 



q^ v = arf" + 



(17) 



This form of the metric requires that the tensor <&'" / must 
satisfy the condition 



71$, 



(18) 



Such feature allows us to write the inverse metric simi- 
larly to the binomial form of the metric, avoiding diffi- 
culties with infinite series 3 . Two remarkable examples 
of this property are the scalar field (in which <fr^ u = 
dpifidvcj)) and the electromagnetic field (in which $ Ml , = 
F a F ) 



A. Special case 

In this section, we limit our analysis to the simplest 
form by setting = u v . In this case, the coefficients 
of the covariant and contravariant forms arc related by 



A=i, 
a 



B = — 



where we set u IJi u L ,r) t11 ' = 1 and write the metric in the 
form 

The associated covariant derivative is defined by 



II f 



where the corresponding Christoffel symbol is con- 
structed using q^ v . The description of an accelerated 
curve 4 in the flat space-time as a geodesies in the metric 
Qfiv is possible if the following condition is satisfied 



0, 



(19) 



where we have used the metric q^ v to write u p 
(a + /3)w M . Therefore, 



3 This is the case of general relativity as a field theory formulation. 
The exact expression for metric tensor is set 



IV. BINOMIAL METRICS 

In recent years an intense activity concerning features 
of Riemannian geometries similar to those described by 
Gordon approach has been done @. In particular, those 
that allows a binomial form for both the metric and its 
inverse. That is, its covariant and the corresponding con- 
travariant expressions are 



(16) 



A consequence is that its inverse, the covariant tensor g M „ is an 
infinite series: 



This formulation was introduced by Feynman, Gupta and others 

(cf. nj). 

4 Note that we are dealing with a collection of paths Y that is 
usually called a congruence of curves. It is understood that each 
element of this collection concern particles that have the same 
characteristics. For instance, if the acceleration is due to an elec- 
tromagnetic field, all particles of T must have the same relation 
between its charge and mass, to wit a bunch of electrons. 
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(tw,-r%,« e ) «" = o. (20) 

In order to preserve the norm u^u^ along the curve we 
assume fi.^u^ = 0, without loss of generality (it corre- 
sponds to a simple re-parametrization along the curves). 
Once the acceleration in the background is defined by 
a fj, = u fi,u u "y the condition of geodetic motion in the 
q^ v -geometry takes the form 

<ty = r 1^u e u u . (21) 
The Christoffcl symbol reduces to 

f%,u e u"= ^-u a u v q av ^. (22) 

Using the expression of q a p in Eq. (|22p and substituting 
the result into the condition ([2~Tj) . it follows that 

°" 2 (a + 0) • 

It means that the acceleration vector a M must be a gra- 
dient of a function VP, i.e., 



only on the angle k a v a between the wave vector k a and 
the dielectric four- vector v a . This is achieved by gen- 
eralizations of the last section. Let us list some examples: 

Case A : the metric q^ v is given by 
rf v = aaf u +phPW, 

and its inverse is 

a -l n P - k k 

a a{a + pN) 

Once the wave vector fc M is a gradient of a given hyper- 
surface £, then we have 

k^v] = 0. 

Substituting this result in Eq. J3|), it follows that k^ must 
satisfy 

= ° ^(g) = const., 

where we define JV( g ) = q^k^k^. That is, in order to 
follow a geodesic motion in q^, k^ must have constant 
norm in the metric q^ v . The explicit expression for this 
constraint is 



a„ = d^. (23) 

Thus, the expression of the coefficients a and (3 of the 
metric q^ v are given in terms of the potential \& of the 
acceleration by 

a + /3 = e- 2 *. (24) 

This simple example gives a very useful formula, which 
exhibits the connection between geometrical and mechan- 
ical quantities. Later, in this paper, we shall analyze the 
case in which VP represents the gravitational potential. 

B. Polynomial Metrics 

Gordon approach depends explicitly on the velocity v a 
of the dielectric. Nevertheless such form of introducing 
an effective metric is not unique. Indeed it seems rea- 
sonable to present another metric q^ v that describes the 
same results obtained by Gordon and besides reduces 
the dependence on the four- velocity v a . From practi- 
cal reasons, it might be useful to weaken this constraint 
of Gordon approach, once it is easier to determine the 
shape of the electromagnetic wave packet in the labo- 
ratory than constructing nonlinear dielectric media with 
arbitrary tensorial parameters e a p and \i a p — despite of 
the great advances in this research area recently [3, EH ■ 

Let us now show that exists a class of geometries 
which play the same role as Gordon metric depending 



N (q) = (a + 0N)N = 1. (25) 

Note that this approach transforms the non- 
normalized wave vector fc M in Minkowski background in 
a normalized time- like vector in q^. It does not vio- 
late Lorentz invariance, because everything happens in- 
side the dielectric. We note that it is not possible to 
fix TV equal to zero, otherwise the metric is ill-defined. 
Therefore, fc M is not a null-like vector in the Q-metric. 
Another feature is that the magnitude of the dielectric 
four-vector 

q^v^Vy = a + P{k.v) 2 , 

is not necessarily positive definite allowing observers with 
velocity great than speed of light inside the medium 5 . For 
instance, if we set 

a + (3{k.v) 2 = 0, 
then, using Eq. (|15[) . wc obtain 

= pipe - l)a 

1 N 



In the laboratory, the angle between these two vectors is easier 
to manipulate than the dielectric velocity field only. We expect 
that this fact could be of reasonable utility in the research of 
analog models. 
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Substituting this result in Eq. (|25|) . yields 

1 

Therefore, the metric q^ v with these values of a and 
/3 produces the following outcome: the wave vector 
becomes a normalized and time-like vector, while the 
dielectric velocity v 11 , which was a time-like vector in the 
Minkowski background, becomes a null geodesic in q^ u . 
Therefore, the causal structure is no more determined 
by fc M . 

Remark that the metric q^ 11 presented in the precedent 
sections is not unique. We can enlarge the set of metrics 
that have the same properties showed above adding 
other terms to <?' il/ provided the condition (fTSj) is valid. 
To exemplify these cases we consider: 

Case B : the polynomial metric is given by 

fn^ = rf v +pk> 1 k u + 8a> 1 a v . 

It is straightforward to show that its inverse has an extra 
term 



N {m) = (I + (3N)N + -N 2 = 0. 

Note that this approach permits a null geodesic 
motion for the wave vector k^. This is the simplest case 
in which we regain the main Gordon result (fc M as a null 
geodesic). The sign of the norm of w' 1 is undetermined 
and may be chosen equal to zero, as we saw in the 
previous case. 

Case C : the most general case involving first order 
derivatives of k^ occurs when the metric is expressed in 
the form 6 

ft* = a t^" + Jfc" A" + s a v a u + Xa^k^ 
and its inverse is 

= -^Hnv + B fc M k v + Aafj,a v + Ka^k u y 

The covariant metric coefficients are given by 
j3(a - Sa 2 ) + Xa 2 



rhtiv = Viiu + B kfj, k v + Aa^a u + Aa^k v ), 

where ( ) means symmetrization. The coefficients of the 
inverse metric are 



B = — 



/3{1 - Sa 2 ) 
X : 



A 



and 



A 



6(1 + PN) 
X 



P8N 
2X ' 



Here, we defined a 2 = —a a a a , N = Nnk^ and 



X = 1 - 6a 2 + f3N - (35 I — + Na 2 



The appearance of an extra term also happens with 
the inverse metric when we consider instead of a^a v a 
term of the form a^k v \ In both cases an extra term 
is necessary breaking the polynomial symmetry between 
the metric and its inverse. Nevertheless we will present 
the calculations for this case focusing only on the metric 
containing the term a^a v and indicating that the results 
are very similar when the other term is considered sepa- 
rately. 

The geodetic motion condition for the wave vector 
leads to 



A = - 



S(a + (3N) - A 2 TV 



and 

where 
Z = a 



A 



X{2a + NX) - 2/3SN 
2Z ' 



N 2 

a8a + a/3N + aNX - {/35 - A z ) ( — + Na' 



Once it involves more degrees of freedom, we can re- 
gain all outcomes presented before, but with different 
algebraic relations. In particular, the magnitude of the 
wave vector in h^ v is set 



N (n) = {a + [3N + XN)N + -N 2 . 

Remark that the metric and its inverse have the same 
number of polynomial terms as required from the begin- 
ning. It did not happen in the case B where an extra term 
was necessary in the inverse metric expression. Following 
this reasoning, in the next section we shall use only the 
cases A and C which satisfy the conditions (fT6f and (|17|) . 
Moreover, as an example, we will set the potential \l/ as 
being the Newtonian potential. 



6 If we use higher derivatives of fc M greater than that which appears 
in the dynamics, then the metric tensor q^y is ill-defined. 
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C. Application: identifying 9 with the 
gravitational potential 

In the weak field limit the description of Newton's 
gravity can be formulated in terms of a geometric repre- 
sentation of the gravitational field. It can be done mak- 
ing use of effective potentials, which correspond to the 
well-known parameterized post-Newtonian approxima- 
tion (PPN) [lf|. In this section, we compare some PPN 
results — particularly, those concerning general relativity 
(GR) predictions, which improve the Newtonian theory 
of gravity in the solar system — with some q— metric given 
by the free-falling (geodesic) condition. We stress that 
q^ v just mimic the geodesic motion characteristic of the 
solutions of GR. Note that there is no dynamical equa- 
tion for q^y and we are not proposing such theory. It is 
purely a kincmatical analogy. 

This section shows that the analysis through geodesic 
paths are much more general than GR, once it is noth- 
ing but a choice of the metric. According to Poincare's 
ideas upon geometrical descriptions of the world fl7| j : 
non-Euclidean geometry is as legitimate as our ordinary 
Euclidean space; the enunciation of the Physics in this 
modified geometry would become more complicated, but it 
still would be possible. In other words, it is possible to 
choose the metric of the space-time, such that an accel- 
erated motion in a given geometry can be described as 
geodesic in another one. 

For convenience, consider Minkowski space-time in 
spherical coordinates 

ds 2 = dt 2 - dr 2 - r 2 dQ 2 , 

and an observer field u M = (1, f(r), 0, 0), which is a gra- 
dient Un = c^E, where £ = t ± F(r) 7 . This vector has a 
non-null acceleration given by 



a il =u ll , v u v = -N^ = (0,-ff, 0,0), 

where prime ' means derivative with respect to radial 
coordinate r and N = u^u^rj^ . 

We choose the scalar function 'I', which characterizes 
the acceleration a M = d^, as identified with the New- 
tonian potential. Then, it follows a relation between N 
and "J given by 

N = 1 + 2* = 1 - — , 

r 

where rn = 2M and M is the mass source of the 
gravitational field (G = c = 1). To go further in the cal- 
culations, we basically split the analysis in two distinct 



This situation can perfectly be adapted to describe wave vectors 
in material media. 



cases. One of them gets immediately a wrong linear 
regime compared to GR, while the other case, which is 
separated in two subcases, can give the expected weak 
field regime. 

Case I : consider the g-metric as follows 

= if +^"/. (26) 

The inverse metric is 

The condition for to follow a geodesic motion in 
this metric is provided by 

u„ ]]u u u = -JV( 9)iM = 0, (28) 

where ii^ = q^ v u u = (1 + (3N)u^ and the magnitude of 
in q^ v -metric is 

N (q) = = (1 + f3N)N, 

which is imposed by Eq. ([27]) to be constant different 
from zero. For convenience, we set N( q ) = 1. 

A power law expansion in terms of e ~ rjj jr of metric 
(f2"T|). which corresponds to the weak field limit, gives the 
following expressions 

goo « 1-^(1 + ^) +G(6 3 ), 

m « -m 3/2 (i + ^)+<^ 7/2 ), (29) 

qn « -l-(?±) 2 (l + UL) + 0{e% 

We see that this metric does not correspond to linearized 
Schwarzschild solution (even in Painlcvc-Gullstrand 
coordinates due to the power 3/2 instead of 1/2 in 
<joi). This metric is similar to some post-Newtonian 
approximation if we consider a rectilinear moving source 
for the gravitational field. The angular components of 
the metric arc identical to Minkowski ones. 

Case II : let us consider another geometry n^ v given by 

Ufa, = rjfiu + B Up u v + A a M a v + A a (p u„) . (30) 
The metric components are explicitly written as 

h 00 = 1 + B, 

noi = f(B — A/'), ( 31 ) 
nn = -l + / 2 (B + A/' 2 -2A/'). 
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The other spatial components are identical to Minkowski 
metric in spherical coordinates. The condition which led 
to follow a geodesic motion in this metric is imposed 
on its magnitude 



N, 



{n) 



where Nua = 



N — A(a 2 N + N 2 /A) 



const., (32) 



i^u v '. We also define 



Z 



NA 



-a 2 A+NB+a 2 NA 2 -BA a 2 N 



N 2 



In this case, we can set either iV( n ) 
us analyze both cases separately: 



or 2V (t[ ) = 1. Let 



Case Il.a: if N, 



that 



( n ) = 0, the assumption (|32l) implies 

N 



A = 



a 2 N + N 2 /4 

Once is null-like in n— metric, we expect to regain 
the metric component of the PPN approximation respon- 
sible to describe correctly the light propagation deflec- 
tion, which is given by the fin component of the metric. 
So, using the considerations above, we set 



noi 
nn 



r 



r 



0(e 2 ), 



(33) 



0(e 3 ). 



Remark that these calculations was basically done in 
order to illustrate some analogies between these geome- 
tries and those studied in analog models of gravity. For 
this reason, the first order approximation is enough to 
show their strong correlation. If someone takes this ap- 
proach looking for a perfect kinematical analogy between 
Schwarzschild metric and some n^, then the introduc- 
tion of a nonlinear scalar potential "J in Newtonian equa- 
tion of acceleration becomes necessary and, therefore, it 
can reproduce exactly the Schwarzschild geodesies, for 
instance. Notwithstanding, we will not enter into the 
details of this generalization because it involves a com- 
plicated question about the physical meaning of the dy- 
namics of such non-linear potential whereas in this paper 
we want to discuss only kinematical properties of the par- 
ticle trajectory. 



V. CONCLUSION 

In this paper, we presented an extension of Gordon 
metric constructing a larger class of geometries which 
describes accelerated motions in Minkowski space-time 
as geodesies in an effective geometry q^ u . This effec- 
tive metric depends only on the background metric and 
on the velocity vector of the accelerated body. In par- 
ticular, we analyzed accelerated paths of light inside a 
moving dielectric and constructed a class of geometries 
with peculiar properties in comparison to Gordon's ap- 
proach, but with similar kinematical effects. Ultimately, 
we gather this new class of geometries, that we call Q- 
metrics, with the effective geometries of nonlinear elcc- 
tromagnetism, producing a collection of possible metric 
structures of the space-time, which has several applica- 
tions in the theory of analog models of gravity. We will 
come back to this in the future. 



Note that a strange linear term appears in noo due to our 
assumptions. If we try to avoid this term making some 
coordinate transformation, then the asymptotically flat 
regime is lost by other metric components, in such way 
that the trouble persists. 

Case II. b : in the case of Nr n \ = 1, we can correctly 
reproduce the linear approximation of Schwarzschild so- 
lution. That is, 

nQ0 « !_Zi + 2(^) 2 + 0(e3), 

noi « 0, (34) 

riii « -l-^ + 0(e 3 ). 

Note that noo differs from PPN results already in second 
order approximation in e (whose the expected value in 
e 2 is 1/2 instead of 2). Therefore, it will surely produce 
some discrepancy in high order terms of the expansion. 



APPENDIX: Hadamard's Method for 
Discontinuities 

We analyze the discontinuities of the electromagnetic 
field according to the standard Hadamard method and 
obtain the dispersion relation for the wave vector fc M . Let 
E be a surface of discontinuity of the field . The dis- 
continuity of an arbitrary function / is given by: 

[/(x)] E = lim (/(z + e)-/(z-e)). (35) 

The field and its first derivative d v A^ are continuous 
across E, while the second derivatives present a disconti- 
nuity: 

= 0, (36) 
MJs = 0, (37) 
[dadpAfj]^ = k a kpCn(x), (38) 

where = c^E is the propagation vector and ^(x) is 
the amplitude of the discontinuity. Substituting these 
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discontinuity properties in the equation of motion 

rf v F av , a = r) av A M;a = 0, 
it follows that: 

k a k p r] af! = 0. 

This means that the discontinuities of the electromag- 
netic field propagate as null geodesies in the Minkowski 
metric -q^. 
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